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Abstract 

We analyze asymptotically a differential-difference equation, that arises in a Markov- 
modulated ffuid model. We use singular perturbation methods to analyze the problem 
with appropriate scalings of the two state variables. In particular, the ray method and 
asymptotic matching are used. 
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1 Introduction 

Fluid models have received much recent attention in the literature. They have been used 
to model statistical multiplexers in ATM (asynchronous transfer mode) networks [TT] . [12J, 
[T5] . packet speech multiplexers [16], buffer storage in manufacturing models [19], buffer 
memory in store-and-forward systems [5j and high-speed digital communication networks 
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[9]. In these models the queue length is considered a continuous (or "fluid") process, rather 
than a discrete random process that measures the number of customers. These models tend 
to be somewhat easier to analyze, as they allow for less randomness than more traditional 
queueing models. 

The following is a description of a fairly general fluid model, of which many variants 
and special cases have been considered. Let X(t) denote the amount of fluid at time t in 
the buffer. Furthermore, let Z(t) be a continuous-time Markov process. The content of the 
buffer X(t) is regulated (or driven) by Z(t) in such a way that the net input rate into the 
buffer (i.e., the rate of change of its content) is r\ [Z(t)}. The function rj (•) is called the drift 
function. When the buffer capacity is infinite, the dynamics of X(t) are given by 

dX = ( rj[Z(t)}, X(t)>0 

dt \ max {tj [#(*)], 0}, X(t)=0 " 1 1 

The condition at X(t) = ensures that the process X(t) does not become negative. When 
the buffer capacity is finite, say B, the dynamics are given by 



, f v[Z(t)\, 0<X(t)<B 

max {77 [£(*)], 0}, X{t) = 
min {rj [Z{t)\ , 0} , X(t) = B 



dt 



The condition at X(t) = B prevents the buffer content from exceeding B. 

In many applications, the process Z(t) evolves as a finite or infinite state birth-death 
process. The description of the motion of Z(t) is as follows: the process sojourns in a given 
state k for a random length of time, whose distribution is exponential with parameter Xk+i^k- 
When leaving state k, the process enters either state k + 1 or state k — 1 with probabilities 

k — > k + 1 w.p. - , k G Jx 

Afc + fJ>k 

k -> k - 1 w.p. ^ k , k e M. 
Afc + /ifc 

The motion is analogous to that of a random walk, except that transitions occur at 
random rather than fixed times. The parameters A^ and /j& are called, respectively, the 
birth and death rates. We shall assume that the birth and death rates are positive with 
the exception of the death rate /xo in the lowest state and (in case of a finite space Af = 
{0, 1, ... , N}) the birth rate Aat in the highest state, which are equal to zero. Also, it will 
be convenient to interpret A^ and as zero if k ^ N '. 

If the buffer has emptied at time t, it remains empty as long as the drift is negative. We 
let rj [Z (t)] = rfe, given that Z(t) is in state k. We shall assume throughout that 7^ for 
all states. We shall also assume that > for at least one k G A/", since otherwise, in the 
steady state, the buffer is always empty. 

We let 

k— 1 
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where an empty product should be interpreted as unity. The stationary probabilities p k of 
the birth-death process can then be represented as 



Pk = , k e Af. 



*3 



When the capacity of the buffer is infinitely large, in order that a stationary distribution 
for X(t) exists, the mean drift Pk r k snou ^ ^ e negative or, equivalently, the following 

stability condition should be satisfied 

$> fc r fe < 0. (2) 

fceA/" 

We let 

Af + = {k e Af \ r k > 0} , Af ~ = {k e Af \ r k < 0}, 
N + = \Af + \, N-=\U~\ 

and since we assume that the drift in each state is nonzero, we have Af + U Af ~ = Af. 
Setting 

P k (t, x) = Pr [X(t) < x, Z(t) — k] ; t, x > 0, k E Af, 
the Kolmogorov forward equations for the Markov process [X(t), Z(t)] are given by 

dPk dPk 

~~dt r k ~dx~ = ^ k ~ 1 ^ >k ~ 1 Vk+iPk+i — (Afc + /i/t) Pk, k E Af. 
For the stationary distribution 

F k {x) = lim 4 _ 00 P fc (t,x) 

we have 

r k F' k = \k-\Fk-\ + fik+iPk+i - (A fc + /i fc ) F k , k e Af. (3) 

Since the buffer content is increasing whenever the drift is positive, the solution to (j3J) must 
satisfy the boundary conditions 

F k (0) = 0, keAf + . (4) 

This means that there is no probability mass at x = if the drift takes you away from the 
boundary. Also, we must have 

F fc (oo) =p k , ke Af, (5) 

as this is the marginal distribution of the regulating process Z(t). In the finite capacity case 
we have the additional boundary condition 

Fk{B)=p k , keAf-. (6) 
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This means there is no probability mass at x = B if the drift moves the process from this 
boundary. The values of Fk(0) for k E Af ~ , and of F^{B) for k G M + , are not a priori 
known. The "half boundary conditions (j4j) and (El) make these problems difficult. 

The purpose of this paper is to continue our asymptotic analysis of fluid models using 
the ray method [6], which we successfully applied in [I] to the model first studied by Anick, 
Mitra and Sondhi in [2]. 

The paper is organized as follows. In Section 2 we state the basic equations. In Sections 
3-7 we analyze these in various ranges of the state space f|T2|) . In Section 8 we study the 
marginal distribution. Finally, in Section 9 we summarize and interpret the results. 

2 Problem statement 

Let Af = {0, 1, 2, . . .}, and the parameters and pk be constant, 

. f a. 1 < k A 

A fc = A, ^ = | Q) k = Q , P = -< L 

The drift is taken as r k = k — c, where c represents the output rate of the buffer. We assume 
c to be a positive non-integer number. This model corresponds to a fluid model modulated 
by the standard M/M/l queue. 

The forward Kolmogorov equations for Fk(x) are then 

(k-c)Fi(x) = \F k . 1 {x)+ t jiF hhl (x)-{X + n)F k (x), 0<k (7) 

fxF Q (x) = AF_ x (x). (8) 



with boundary conditions 
and limiting distribution 



F fc (0) =0, LcJ + 1 < k. (9) 



F fc (oo) = (l-p)p fe , 0<k. (10) 
Here (jHJ) defines F_i(x) and this condition is equivalent to 

-ci^(x) = /^(x) - AF (x). 

Since the buffer capacity is infinite, we need the stability condition 

< c, or p < 1 —. (11) 



1-p ' c+1 

A related model, with r Q = p < and r k = p > was studied in [18], [3] and [1] where a 
spectral representation of the solution was obtained. The same model was analyzed in [10] 
using continued fractions. The general case, with arbitrary r^, p k and A&, was studied in 
[7] and [IT] using a family of orthogonal polynomials. The fluid queue driven by a general 
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Markovian process was analyzed in [13]. A numerical method was presented in [§]. The full 
transient solution was considered in [13]. 

We shall analyze this model directly by using the differential-difference equation (J7J satis- 
fied by Fk(x). After appropriate scalings of k and x, we analyze this equation asymptotically 
for c — ► oo, using singular perturbation methods. We also carefully treat various boundary 
and corner regions of the state space 



{(x,k):x>0, 0<k} 



(12) 



and indeed we show that their analysis is needed in order to obtain the asymptotic expansions 
away from the boundaries. 



3 The ray expansion 

To analyze the problem (171)- (11 01) for large c we introduce the scaled variables y and z, with 

x = c 2 y, k = cz, z,y = 0(1). 

We define the function G(y, z) and the small parameter e by 

e = c -1 , F k (x) = G [xe 2 , ke) = G(y, z) 

and note that F k ±i(x) = G(y, z ± e). 

Then (J7J becomes the following equation for G(y, z) 



dG 

e(z - 1 )-q-(Vi z ) = AG (^' z ~ £ ) + z + e)-(X + fi) G(y, z) 



G(0,z) = 0, l<z. 



and ([nD implies that 
Also, from ffTOj) . we have 

F k (oo) = G(oo, z) = (1 - p) exp 



-^ln (p) 



< z. 



(13) 
(14) 

(15) 



To find G(y, z) for e small, we shall use the ray method. Thus, we consider solutions 
which have the asymptotic form 



G(y, z) ~ e u exp 



K(y,z), 



(16) 



where v is a constant that must be determined. Using (flBj) in (fl3l) . with 

-*(3/, z ± e) = ± * 2 + Ul zz e + O (e 2 ) , 



dividing by exp [-^(y, z)\ , and expanding in powers of e we obtain the eikonal equation 
for ^(y, z) 

^(l-e q ) + \(l-e~ q ) + (z-l)p = 0, (17) 
and the transport equation for K(y, z) 

- Ae"«) ^ + (1 - z) ^ + ig + Ae-) X = 0, (18) 

where 

_ dm _d^ 

^ dy' dz 

To solve (jTTJ) and ( fTBl) we use the method of characteristics, which we briefly review below. 
Given the first order partial differential equation 

3 (y,z } ty } p } q) = 0, 

where p = q = we search for a solution ^(y, z). The technique is to solve the 

system of "characteristic equations" given by 



y — a , — Bp, z — $ q 



dy_ 

dt 

P=-$ y -p$y, q = -$z-q$v 

?P = p$p + q$ q 



where we now consider {y, z, ip, P, q] to all be functions of the variables s and t, with ip(s, t) = 
^(y, z). Here t measures how far we are along a particular characteristic curve or ray and s 
indexes them. 

For the eikonal equation ( 1TT1) . the characteristic equations are 

y = z - 1 (19a) 

z = \e~ q - ne q (19b) 

p = (19c) 

q = -p (19d) 

iP = p( z - 1) + q (\ e ~ q - /ie 9 ) . (19e) 

The particular solution is determined by the initial conditions at t = 0. We shall show that 
for this problem two different types of solutions are needed; these correspond to two distinct 
families of rays. 

Setting ^ y \ t=0 = s, ^z\ t=0 = B and solving ffT9cD - fll9dp . yields 

p = s, q = B-st (20) 

so that ty y is constant along a ray. 
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3.1 The rays from (0, 1) 

We now consider the family of rays emanating from the point y — 0, z — 1. Evaluating 
(P2D at t = we get 

// (l — e B ) + A (l — e- B ) = 

so that 

B = or 5 = In (p) . (21) 
From (119bl) and (|20|) . with the initial condition z(s,0) = 1 and using (l2Tj) . we obtain 

2 = i [Ae" B (e si - 1) + fie B (e~ st - l)] + 1. (22) 

From (|19ap . we have 

y{s, 0) = z(s, 0)-l = 

and 

y(s, 0) = z(s,0) = \e- B - fxe B . 

From fl2Tl) we have 

A-//< 0, fi = 



^ (S ' 0) \ /x- A> 0, B = ln(p) • 
Using the initial condition y(s, 0) = and expanding in powers of t, we get 

t 2 

y(s,t) ~ y(s,0)-, 

and in order to have y > for t > (i.e., for the rays to enter the domain [0, oo) x [0, oo)) 
we need to choose 

B = In (p) (23) 

with £> < since p < 1. 

Integrating (II 9 aft and using (T22]) and fT23|) . we conclude that 

y(s, t) = ^[fi (e st - st - 1) + A (l - st - e~ st )] (24) 

z (s, t) = - s [fi (e st - 1) + A (e~ st - 1)] + 1. (25) 

This yields the rays that emanate from (0, 1) in parametric form. Several rays are sketched 
in Figured! 

For t > and each value of s, (|24|) and (|25|) determine a ray in the (y, z) plane, which 
starts from (0, 1) at t — 0. We discuss a particular ray which can be obtained in an 
explicit form. For s = we can eliminate t from (1251) and obtain 

( z ll^ 

y = Y {z) := 2 (~_ A) » * = 0, 1 < (26) 
7 




Figure 1: A sketch of the rays from (0, 1). 



and along this ray, t and z are related by 



t{Y ,z)=T {z) 



z - 1 
\i — A 



, 8 = 0, 1<Z. 



(27) 



For s > 0, we have both y(s,t) and z(s,t) increasing for t > 0. For s < the rays 
reach a maximum value in z at t = T\, where 



T 1 ( S ) = ^-ln(p), s<0 



and we have 



y(s,Ti 



1 



+ 1. 



2vA/i - (A + fi) 



(21 



(29) 



From (I19al) we see that the maximum value in y is achieved at the same time that z = 1, 
and that occurs at t = T 2 with 



T 2 (s) = -ln(p), s<0 

s 



(30) 
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and 

S/( S ,T 2 ) = i[2(A-/x)-(A + /i) ln(p)]. 
Inverting the equations (1211) - (1251) we can write 

s = S(y,z), t = T(y,z) 

and 

2) = -0 [5 (2/, ^) , T(y, z)] , K(y, z) = K[S (y, z) , T(y, z)} . 
We will use this notation in the rest of the article. 

3.2 The functions * and K 

From (j!9ep we have 

V> = /ie st [1 + In (p) - ts] + Xe~ st [1 - In (p) + ts] , 
which we can integrate to get 

Ms, t) = ^e st [2 + In (p) - ts] - -e~ s * [2 - In (p) + ts] (31) 

s s 

+ ^M)-^[2 + ln(p)] + ^[2-ln(p)]. 

Obviously, ip(s, 0) = ipo is a constant, since all rays start at the same point. Setting s = 
in ( 13TT) and using ( 1271) . we obtain 

^(0,t) =Vo + (*-!) I* (/>) 
and therefore, taking the limit as t —>■ 00, we get 

*(oo,z) = Vb + (*-!) to (p)- 
On the other hand, from ( 1151) we have 

^(00, z) = zln (p) 

and we conclude that 

i) = In (p) . 

Solving for e st in (]24l)- (T25T) . we get 

e st = l + ^[z-l + y8 + t(X + fi)], e~ st = l + ±.[ z -l- y8 -t{X + fjL)]. (32) 



Replacing ( 1321) in ( 1311) . we obtain 

i) = 2ys + [In (p) - st] (z - 1) + In (p) . (33) 
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We shall now solve the transport equation (Tl8l) . which we rewrite as 



dy 

Using (El) and (ESI) in (El), we have 



:,-i)* + (^-^«E = i!o«. + ^K. 



dK ldq 



(/ie 9 + Xe~ q ) K. 



dt 2dz 

To solve ( 1351) . we need to compute || as a function of s and t. Use of the chain rule 



'dy 


dy' 




-dt 


dt- 




"1 


0" 




ds 




dy 


dz 




$z 


dz 




a| 


ds 







1 


.dt 


ds. 




-dy 


dz. 





-at 


at- 


1 


dz 


dy' 


dy 


dz 


ds 




sL 


ds 


" J 


dz 


dy 


-dy 


dz. 


L at 


dt J 



dz , n9 
-(*-!)' 



and hence, 



where the Jacobian J(s, i) is defined by 

. ch/ <9z <9y dz 1 
<9i 9s 9s dt s 

Using (1551) we can show after some algebra that 

dq _2y 
dz J ' 

while (1321) gives 

/ze 9 + Ae~ 9 = // + A + (z - l)s. 
Thus, the transport equation (155]) becomes 



1 97^ u r x . „. , 

if^r = -j^ +A+ ^- 1 w 



Using (|23j) and (|25|) in (|37|). we have 

dJ _ 2y d 2 z 
~dt ~ ~7W ~ 

Combining (|39|) and (HOj) . we obtain 

1 <9K 



if at 



whose solution is 



K(s,t) 



2y\jM + \+(z-l)s] 



_ 1 dJ 
~2J ' 
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where K (s) is a function to be determined. 
From (1371) . we have 



J( S ,t) = (/i 2 -A 2 )^ + 0(t 4 ), f-0. (42) 

Since the Jacobian vanishes as t — > 0, the ray expansion eases to be valid near the point 
(0, 1), where a separate analysis is needed. 

So far we have determined the exponent ip(s,t) and the leading amplitude K(s,t) 
except for the function K (s) in (|4ip and the power v in ffTB"]) . In Section 4 we will 
determine them by matching (fl~6l) to a corner layer solution valid in a neighborhood of the 
point (0,1). 

3.3 The rays from infinity 

Denoting the domain in the (y, z) plane by 

£> =[0,oo) x [0,oo), (43) 

we must determine what part of D the rays from infinity fill. The expansion corresponding 
to these rays must satisfy the boundary condition (IT51) . Thus, we have 



p(oo, z) = poo = 0, g(oo, z) = goo = In (p) , (44) 
while fl!9ap - (119bf) yield equations for the rays yoo{t), z^t), 

yoo = Zoo - 1, ioo = /i - A (45) 
or, eliminating t from the system (145]) and writing y^t) = Y^z) we get 

dYn z - 1 



cfe fi — X 

Solving (|45l) subject to the initial condition Y^zq) = y , where y x z = 0, we get 



(46) 



Yoo(z) =Vo+ 2 ^_ x) K z - !) 2 - («b " I) 2 ] • (47) 

From (1461) . it follows that the minimum value in y occurs when z — 1. Hence, for y to be 
positive, we must have 

where 5 / o( 2; ) was defined in fT26|) . Therefore, the rays from infinity fill the region given by 

R = {0<y, 0<z<l}u{Y (z)<y, I < z} . (48) 
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Figure 2: A sketch of the rays from infinity. 



The complementary region R 



c 



R 



c 



{0<y< Y {z), 1 < z} 



(49) 



is a shadow of the rays from infinity. In R c , G is given by (fl6l) as only the rays from (0, 1) 
are present (see Figure [2]). In the region R, both the rays coming from (0, 1) and the rays 
coming from infinity must be taken into account. We add (fl6l) and (fT5l) to represent G in 
the asymptotic form 



G(y,z) ~ (1 -p)exp 



—z In (p) 



+ e exp 



-$(y,z) 

e 



K(y,z), {y,z)eR. (50) 



We can show that z In (p) > ^f(y, z) in the interior of R, so that G(y, z) ~ G(oo, z). However, 
in R we can write fl50|) as G(y, 2) — G(oo, z) ~ e^exp [^(y, 2;)] 2). 



4 The corner layer at (0, 1) 

We determine the constant v in (TX6T) and the function K (s) in (|4ip by considering carefully 
the region where the rays from (0, 1) enter the domain D, and using asymptotic matching. 
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We introduce the stretched variable I, the function Fj and the parameter a defined by 

F k (x) — G x (x, k — c + a) — 

l = k — c + a, — oo < / < oo (51) 
a = c — [_c| , < a < 1. 

Note that a is the fractional part of c and I takes on integer values. Use of (!5~T!) in (J7j) yields 
the equation 



(I -a) 



dFl 



(i) 



fiF^\ + \F£{ — (X + fj) F^> , x > 0, Z G Z. 



(i) 



Also, (jHD gives the boundary condition 

F/ 1} (0) = 0, 1>1 
and (TT5l) implies that Ffc(oo) = Fj (oo) with 

F z (1) (oo) = (l-p)p ; - a exp 



In (p) 



(52) 



(53) 



(54) 



For a fixed / and x — ► oo, we approach the interior of R, where (1501) applies. Thus, (|54p is 
the asymptotic matching condition between the corner layer and the solution in R. We shall 
examine the matching to R c later. 

Since the problem (152]) - (154"]) is of interest by itself, we solve it in a slightly more general 
setting in the next theorem. 



Theorem 1 Suppose that the function <&i(x) satisfies the equation 



dx 



A® 



i+i 



C$1-! -B$i, x>0, lez, 



with 

and the boundary conditions 
where 



B > 2VAC = (3 > 

$/(0) = 0, 1>1, 
$z(oo) = ^r 1 , leZ, 

Ar 2 - Br + C = 0. 



Then, has the integral representation 



C\ 3 A 1 



.4 



B2wi 



exp ( x9 + 



Br 



x -r 1 - a + 



B 



B-A 




P\ (B- AB 



(55) 

(56) 

(57) 
(58) 

(59) 
(60) 



de 
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and the spectral representation 



, ^ / fC\~ 2 /A <A (j + l-a) 3 , \ 

V / i=0 

x exp + ^) g) , 

where Br is a vertical contour in the complex 8— plane on which Re(6 l ) > 0, T (•) denotes 
the Gamma function, J v {-) is the Bessel function of the first kind, 



A = ^B 2 -(3 2 = VB 2 - AAC 

and 



B 

j + l-a' 



(>;, ' ~ ; - j = 0,1,.... 



Proof. Equation fl55|) admits the separable solutions 

$ l {x) = e 9x h l {9) (62) 
if hi(6) satisfies the difference equation 

Ah l+1 + Ch l „ 1 = [(l-a)6 + B] /i,. 

Setting 

*i(*) = (l) 2 ^)' 

we see that 

+ = I [(/ - a)0 + 5] (63) 



The only solutions to (1631) which have acceptable behavior as / — ► 00 are of the form 

where J v (-) is the Bessel function. If ( 1621) is not to grow as x — > 00, we need 9 < 0. But 
except when v is an integer, the Bessel function </„(•) is complex for negative argument. 
Therefore, we need 

B 

-a + - = -1,-2,... 

or 

6 3 = -—^ <0, j>0. (64) 

j + l-a 
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It follows that the general solution to (155]) takes the form 

QAx) = $/(oo) + 



i>o 



or 



j>0 



exp 



B 



j + 1 - a 



x Ji. 
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(j + 1 - «) 



(65) 



where the coefficients aj in the above (spectral) representation remain to be determined. 
Taking the Laplace transform 



of (1631) we obtain 



*,(0)=/<„/i+ 



7 £ 



j>0 j+l-a 



(66) 



Thus, the only singularities of $z(#) are simple poles at 9 = and 9 = 9j, j > 0. It is 
well known that the Gamma function T(z) has simple poles at z = 0, —1, —2, . . .. Hence, 
we shall represent $/(#) as 



(67) 



where /(6 1 ) is chosen such that 



5 



1 — a 1 J, 



/(*) 



is analytic for Re(6 l ) > — Taking the Laplace transform in (1531) . we get the equation 

(/ - a)9$i(9) = A$ l+1 (9) + C$,_i(0) - B$i(0), I > 1, 
which is satisfied by ( 167]) . By the inversion formula for the Laplace transform, we have 



cy i 

A/ 2^i 



e xe ^r I r. + l _ a JJ, 



1_ 



V 9 



'l-oc+% 



t V 9 



f(&)d9, 



Br 



where Br is a vertical contour in the complex 9— plane on which Re(0) > 0. 
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Since the residue of $i(9) at 9 = corresponds to <&j(oo), we must have 



C 



B 



e 



as 9 — - > 0. Using the asymptotic formulas (z — > oo) 



T (a + 62) ~ v^e" 62 (bz) a+bz ^ , 6 > 



and 



Ja+bz \CZ) 



in ([5711 . we see that 



r exp 1 z 



'Vb 2 — 



b-Vb 



2 _ n 2 



a+bz 



, b > c> 0, 



r 



5 
I" 



1 — a ) J, 



'l-a+ 



f V 9 



or, using (159|) . 
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Therefore, we write 



~ e 9 



B - AY fB - AB 







(3 9 



f B - AB\^~ a IB 



(3 9 



f{0) = ii oo J-exp [T(9)]f(9), 



where 



T(0) 



B-A (B 



9 



9 



-a In 



B-AB 

(3 J 



and f{9) is entire, with /(0) = 1. 

By combining the preceding results, we have 



A' e ^°- 



1^ fC\~ 2 1 f r0 l^ (B 

Br 

The boundary condition (137)1 implies that 



exp [T(0)]/(0)d0. 



lim 



and using the asymptotic formula 

zx 1 ' 1 



2) T(u+iy 



0, Z > 1 
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in (!72|) . we have 



1„ (B 



9 \ 9 



1 - a }J l 



2(3 



exp [T{9)\ 



i+i 



r(l-a) f/3\ l (rB\ a 



9 J T(l-a + l)\2j \C 



, 9 — > oo. 



Setting Z = 1, we get 

f{9) =o(9), 9^oo 

and Liouville's theorem forces f{9) to be a constant. Since /(0) = 1, we have f{9) = 1. 
Thus, ( 1721) becomes 



xT 



'A /C 



B \A 2vri 



5 



27Ti/ 


xB 1 




9 


Br 




Jl-a+f 


(§) 



(73) 



d9 



and we obtain (1501) . 

The coefficients a, in the spectral expansion (1651) are determined from (1731 by applying 
the residue theorem. Noting that 



Res 



1_ (B 



9 \ 9 



+ 1 - a , = 0,- 



'A (C 



we obtain 



Using (1641) in ( J7TT) . we get 



(j + l- a)j! 



A,J b (§- ) exp [T(^)] , j>0. (71) 



exp [T(0j)] = exp 



B-A\ (B-A 



9; 



ft 



3+1 



(a-j-1) 



and (IBTj) follows. ■ 

For the problem (|52l) -( l54l) . we have 



A = fjt, B = X + p, C = X, r = p, /ioo = (1 -p)p c - a , A = p-X (75) 
and therefore 



F/ 1) (x) = (l-p) 



/i + A 27ri 



xT 



X + p 



+ 1 - a) J, 



Br 







(76) 



exp [A(9)) 



d9. 
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with 



Also, 



. 2A /A + p \ { \ + n 



F/ 1) (x) = (l-p)p c -^ 



H - A v-^ (j + 1 - a) 



E 



(77) 



xp^ 1 exp ( xtij + J,_j_i \^J^ 



where 



X); 



A + p 



J >0. 



j + 1 - a ' 

This completes the determination of the spectral and integral representations of F[ (x) 
and hence the leading term for Fk(x) in the corner region. 

4.1 Matching the corner and R c regions 

In this section we shall determine the function Kq(s) in (j4ip and the power v in f fl6|) . We 
begin with the following result. 

Theorem 2 With the same hypothesis and notation as in TheoremUl let Q be defined by 

2Ax 



n 



0(1). 



(7f 



Then, 



25 



(1 — fl) , x — > oo, Z — > oo, f2 — ► 1. 



3ttA(Z - a) 
Proof. We set 

# = e e, 77 = (Z-a)0 + B = (z-l)6 + fl, 77,0 = 0(1), 77,6 >0 
and using (I69I) . we obtain 



-exp I —y/v 2 -P 2 



i 



»7 



n 

_-<-> 



ft 



eQ 



2vrp(p) 6XP I eQ 



1 



p (,) + ,ln(^) 
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with 

Use of (1681) gives 



r (l + 1 - a 
and from flTTj) we have 

exp [T(e©)] = exp 

Therefore, 



27rex P ^ 



B-A B , 



B-AB 



(3 9 



-) 



B-AB 







5 



J *(^© ) r (^ + 1 -«)^P [T( £ 0)] 



B f B — A 



e6 
p(r?) + ryln 



Using l j80l in ( |73l) yields, in terms of z and Q 

i /■ 1 



y-p(v) \ 
P ) 



A -Bin 



B-A 



A/ 27riJ V-B J^nj 

Br' 



exp 



(z - 1)0(77) 



where 



. . (rj-B)Q 1 



2A 



77 - B 



Pi 7 ]) + Tiln 



77 ~P(^) \ 
/ 



A - Bin 



B-A 



(79) 



50) 



51) 



52) 



and Br' is a vertical contour in the complex plane with Ke(r]) > B. For e — > 0, with f2 fixed, 
we can evaluate (1BT1) by the saddle point method to get 



Hx) ~ /W"VA ( - 



z-l 



2tt (z-1)?7*-B 



exp 



-{z-l)g(v m ) 



53) 



where the saddle point 77* (Q) is defined by #'(77*) = 0. Note that 77* (fl) > B for fi < 1, i.e., 
the saddle point 77* (O) lies to the right of the pole at 77 = B and the integrand is analytic 
for Re (77) > B. 

Taking the derivative of ( 1521 . we find that 



s'fo) = — + 



2A - B) 



A + Bin 



B-A 

77 - p(?;) 



-p(rj) 



54) 
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and we observe that g'{ff) = if r]* = B and Q = 1, which implies that rj* (1) = B. To 
determine rf for Q ~ 1, we use and an expansion of the form 

rj* (fi) ~ 5 + fll (O - 1) + a 2 (fi - l) 2 + a 3 (fi - l) 3 + • ■ • . (85) 

Using ( |85i) in ( 18^1) and expanding the latter in powers of Q — 1, we find that 

3A 2 27A 2 
a 1 = -— , a 2 = -^(A 2 -35 2 ) 



and 



5-A\ 3A _ _ 2 £ 



^)~ln^— j-^(^-ir, 3A3 
1 25 ,„ _! 1 1 



(n-i)- 



?7* -5 3A 2V 7 ' v /^yy vS' 

from which we conclude that 



•'W—^'3^o (n " irl - (86) 



Using (J7o| in (job]) , we get, for x, Z — > oo with — > 1 

z In (p) 



~ -(1 -p)exp 



2(jj + \)e 



37r(/x — A)(z — 1 



(n-i)- 1 . (87) 



This must agree with the behavior of the ray expansion in R c as (y,z) (0,1). We next 
evaluate K and ip in (fTSj) near the corner (0, 1). From f[2"5"j) . we have 

, 2 ~ 1 (2:- l) 2 (/i + A) . 

r ~ n — s, s — > . (88) 

/i-A 2(/i-A) 3 V 7 



Using flggp in (E3p. we obtain 



or, using ([78 



* 2( M -A) 3( M -A) 3 ' 



30.-A)'n-i _ L 



2 (// + A) z - 1 
We expand (l31l) for small t 

V>M) ~ In (p) [1 + (// - A)t] , t->0 
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which taking ([88]) into account gives 

>%,z) ~ln(p) + (z-l)ln(p), z 

in agreement with the exponential part of ( J87l) . 
From (HB-tpi). we obtain 



K(s,t)~K {s) 



or, using (189]) and (IHHj) in the above, 



(p 2 - A 2 )t 3: 



t -> 0, 



3 (p - A) 2 Q - 1 



3(p - A) 3 



2 + A) z- I J y (p 2 -A 2 )(z-l) 3 ' 
Matching the algebraic factors in ( 1H7I) and ( 1901) yields 
3(p- A) 2 fi- 1 



z -> 1. 



(90) 



e v Kn - 



2 (p + A) z-1 



2 p + A 2 — 1 //i - A 

£ 3(/,-A) 2 fi-iV^r (1 - p) 



which implies that 



#o(s) 



p — A 1 — p 



(91) 



2vr s 

and v — |. This completes the determination of the asymptotic solution corresponding to 
rays from the point (0, 1). To summarize, we have established the following. 



Result 3 The solution of ( T73l) is asymptotically given by 



G(y, z) ~ v^ ex P 
G(oo, z) — G(y, z) ~ — \fe exp 



K.(y, z) in R 



c 



K(y, z) in R 



(92) 
(93) 



with 



-zln (p) 



G(oo, z) = (1 — p) exp 
# = {0<y, < z < 1} U {Y (z) < y, 1 < z} , Y (z 



(z-1) 



2(p-A)' 

*{y, z) = t) = 2ys + [In (p) - st] (z - 1) + In (p) , 



1< z. 



p — A 1 — p 
27rJ(s,i) s ' 



(94) 
(95) 
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where (y,z) is related to (s,t) by [Efy and &5\) , and J(s,t) was defined in We note 

that s < in R so that the right side of / fPffj) is positive. This gives the leading term for the 
probability 

Pr 



y % 

X(oo) > x = -, Z(oo) = k = - 



that the buffer exceeds x = cy. 

In the corner range where 137]) applies, the leading term is given by (f76|) or ([770 



5 Transition layer 

We shall find a transition layer solution near the curve y = Y (z) defined by f[2"6"j) which 
separates R and R c . On this curve s = 0, and hence f l9li)j) is not valid because K(y, z) is 
infinite there. 

We introduce the new function L k (x) defined by 

F k (x) = F k (oo)L k (x). 

Then ([7]) yields for L k (x) the equation 

(k - c)L' k = /uLfe-! + \L k+1 - (A + fi) L k 

and the boundary condition (fTOl) becomes 

L k (oo) = 1. (96) 

In terms of the variables y = e 2 x, z = ek, the function L^(y, z) = L k (x) satisfies 

e(z - l)^-(y, z) = z-e) + AL«(y, z + e) — (X + fi) L«(y, z) 

dy 

and therefore, as e — >■ 0, 



Introducing the stretched variable g, defined by. 



dz 3 



( z ll^ 

y = Y (z) + ^~eg= ^-Z- > + sfeg (97) 

and the function L^ 2 '(g,z) = L^'(y,z), we obtain for L^ 2 \g, z), to leading order, the 
diffusion equation 

(X a) dL(2) I (A + ^ } (z iV^-O (98) 



22 



To solve (I98I) . we assume that L^ 2 '(g,z) is a function of the similarity variable V = ^y, 
and let £(V) = L^ 2 >(g,z), where r(z) is not yet determined. From (I9"gj) we get 



(A- fi)r(z)r'(z)V£! + 



' 1 (A + /i) ^-l) 2 £" = 



(A-/i) 
and (196]) gives 

fl(oo) = 1. 

We can eliminate z in ( |99l) by choosing r(z) to satisfy the equation 

_(A- A4 )r(^) = -^±4 (^-l) 2 . 

(A -fjt) 



(99) 



(100) 



(101) 



We choose r(l) = 0, which is necessary for matching the transition layer with the corner 
layer solution (JTSj) . and solve fllOip to obtain 



r(z) 



'2 /x + A 
3 (M - A) J 



(102) 



Now flU]) and ffTOOl) become 



£" = _V£', £(oo) = 1 



with the solution 



Z(V) 



1+erfl ^) 



/exp ( — m 2 1 du. 
P l 2 ' 



Thus, the transition layer solution for y — Y (z) = O (y/e) and 1 < z is 

1 



F k (x) ~ (1 - p) / 



1 + erf( ^) 



(103) 



with 



3/-K (*) /3(/x-A) s 



V 2 // • A 

_ y-Y (z) 
r(z)y/e 

to the ray expansion in R c , as given by fl92|) . corresponding to rays emanating from (0, 1 



We can show that if (11031) is expanded as V = v r (J°^ ~^ — oo, the transition layer matches 



23 



6 The boundary layer at z = 



The ray expansion in ( 1931) does not satisfy the boundary conditions in (JHJ). Thus, we re- 
examine the problem on the scale z = 0(e) (k = 0(1)) , with y > 0. We consider solutions 
of ([7]) which have the asymptotic form 



F k (x) - F fc (oo) = F^(y) - F fc (oo) ~ s" s exp 

Using (11041) in and expanding in powers of e gives, to leading order 

= XKfl, + tiKfl + [%(y, 0) - (A + p)) K { *\ 

or, using ( 1201 . 



= XKfl, + nK^ + [S(y, 0) - (A + //)] 



.(3) 



(3) 



and from (JHJ) we get 



pK®=K® 
From ( 1251) we have, along z = 0, 

/^(y) + Ar^Z/) + S(y, 0) - (A + y) = 0, 

where 

f(y)=exp[S(y,0)T(y,0)] 
and therefore we can rewrite (11051) as 



K, 



(3) 



pC\y) + m- 



Since S{y, 0) < and T(y, 0) > 0, we see that < £(y) < 1 for all y. Using (jTQgj) in 
have 

^(2/) " AT%) + A - // - (A + //) In [£(y)] = S 2 (y, 0)y, 
which combined with (11071) gives 

(i - r 1 ) p - a - e) - (p + 1) in (o = a* [(i - r 1 ) p + (i - o] 2 y- 

Solving (I109p subject to (11061) . we obtain 



K?\y) 



m - 1 



e(y) + P k C k (y) 



[i-p^'iy] 

with K(y) to be determined and hence, 



K(y), 



Fl 3 \y)-F k (oo)^e U3 exp 



1 



*M) 



m 1 : e(y) + P k r k (y) 



i-pt l {y) 



(104) 



(105) 
(106) 

(107) 
(108) 

(109) 

we 



(110) 



111) 



K(y). (112) 
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Setting k = z/e, F^\y) - F k (oo) = z) in ffTT^j) and letting e -> 0, we get 



~ ^ 3 exp <j -#(?/, 0) + ^zln 



or, using fl20|) . 



z) ^ e vz exp 



From (JMD, (EZD, and (TTO we have 



1 1 9$ 

0) + -2— (y, 0) 



G(y, z) — G(oo, z) ~ v^ ex P 
as z — > 0, with 



i i a* 



fi — A 



27rJ (y)-S(y,0) 



(113) 



;i - p) , (H4) 



V(y, 0) = 2y<%, 0) + In [£(y)] < 0, J (y) = 2 

o) = (A + /i) - tf{ y ) - \r\y) < o. 

Matching (11131) and (11141) we conclude that 



ay) 



y-l<0, 



^3 = 2^ k{v) 



\i — A 



27TJ ( 2 /)S( 2/ , 0) 

Therefore, for k = 0(1) and y > 0, we have 

Fl?\y)-F k (oo)~ v^exp 



(1-P)- 



-*(2/,0) 



(115) 



where £(?/) is defined implicitly by fll lOj) . 

7 The boundary x = 

For x = and fc < [c\ , the values of F k (0) can be computed from the ray expansion, since 
Ffc(0) -F k (oo) ~ v 7 ^ ; z) ex P [i*(0> *0] is wel1 defined. For x = and fc > [cj -hi, we have 
Fk(0) = by We now examine how this boundary condition is satisfied by considering 
the scale y = 0(e) and z > 1. Note that this part of the boundary is in the region R c . 
From (1251) we have 



^ - 1) s + /x + A + \l [(z - 1) s + /x + A] - 4A/x 



z>l. 



(116) 
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Using (PTBp in ([24]) we get 

2) = + — — {(n + A) In 

1/ 21 

and using (I117p in (11161) 



+ 2\^+0{y\n 2 y), (117) 



r(y,z) 



^-l 



In 



y + 0{y 2 \v?y) 



for |/ — ► + and z > 1. 

Using (flTTD and ffTT8l) in (El and (jSS}, we find that 

V{y, z) ~ 2) = In (p) + (2 - 1) | In 



Ay 



2-1 



(A + n) In 



2-ir 

+ ^-v}y, y 



h2 
0. 



Hence, we shall consider asymptotic solutions of the form 



F k (x) ~ e a exp 



~^(ex,ek) 

£ 



K(u, ek), 



'118) 



(119) 



(120) 



where u = ex, u = 0(1) and cr, K(u,z) are to be determined. Using (I120p in f[T3"j) we get, 
to leading order, 

s dK dK , . 

( z -l) — + u— + K = 0. (121) 



<92 du 



The most general solution to (11211) is 



Hence, 



K(u,z) = -k(E), 
u 



F k (x) ~ G(u,z) = e a exp 



2 - 1 



■it 



: ^(eu, 2) 



:*(S). 



(122) 



(123) 



To find k (S) and o we will match (I123p with the corner layer solution (1761) . 

Recalling that I — a — £ — !- and using the asymptotic formula (1691) we get, as e — ► with 
fixed 



t/ z-l 



2V/^A 



2vr (2 - 1 



'2 - 1 A + fi . , 
cxp <; | — 1 — ) In 



\i\ee 



9 z-l) 



(A + //) 



(124) 
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Using (11241) and writing (!76j) in terms of u = ex and z = 1 + (I — a)e, we have 

z + 1 oT 



X 2^i/{r XP 



Br 



xr 



A + p 



p — A 


exp 


p + A 


— + 








"A 


-/i 



2e 2 
2 — 1 A + p 



ln(p) 
In 



2vr (-2-1) 
pAee 



e 



e 



a In ^/p 



A + p 



(125) 







d0, 



To evaluate (11251) asymptotically as e — > we shall use the saddle point method. We find 
that the integrand has a saddle point at 9 = S, so that 



F«(z)~ £ (l-p) 



I fx - X 
P + A 



exp 



— --)ln(p) 
2e 2 ' v ^ 



1 1 

27TM H 

LiXee 



_ / A + p \ TmS f z — l A + p\ . / y/JIXee \ 



or 



x exp 



x exp 



A - p / 'A + p \ / _ A + it 
— a I In I VP — 



i-p i i r/ a + /i 



p 27TU I 



In (p) us , / Ae e 
— — H In ' 

e e 



S 2 M 



a In 



+ 1 - a 
A + p 



A + p 
x exp < — = — m 



Eu (p + 1) 



2A 
+ — 



Writing (11231) in terms of H, we obtain 
G(w, mS + 1) = e a exp 
x exp 

Matching (T126D with ffT2Tj) . we have 



In (p) , / Ae e 

— — H In ' 



mH 2 



-*(H) 



fc ( 3) = (1 _ ) ,) ( ,i^J_Ir(^±ii + ,-„ 



x exp 



1 + p 27T S 

Mini ^+^ln / ' 



A + p/_ 



(126) 



(127) 
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and cr = l. Therefore, for 1 < z, 



G(u, uE + l) = e(l-p) jl-^ J-ilr f + 1 - a 



1 + p 27T H u 



x exp 



In (p) uE , / Ae 2 e 
— — H In ' 



S 2 M 



a In 



A + /i 



A + p 
x exp < — — — m 



& (p + 1) 



or 



G(u,z)=e{l- p)J^—^- — J' 



In (p) z — 1. 
x exp < — — H In 



1 + p2irz - 1 
Ae 2 e-u 



2A 



(A + p) u 



z-l 



+ 1 - a 



x exp 



e e 
(X + p)u 
z-l 



(z-iy 



+ a In 



In 



LO + i)0-i)j z-i 



+ 



(A + p) u 

z-l 
2Xu 



Note that from ffT2gl) we have G(u, eife) = O (u fe ~L c J) , as u -»• 0, Jfe > [cj + 1. 

8 The marginal distribution 

We will now find the equilibrium probability that the buffer content exceeds x, 



M(x) = Pr [X(oo) > x] = 1 - J^x) 



fc=0 



for various ranges of x. 



(128) 



(129) 



8.1 Approximation for x = 0(1) 

In this region we shall use the spectral representation of the corner layer solution. Using the 
generating function 

oo 

Ji{x)z l = exp 

i— — oo 

in the form 



x f 1 

— z 

2 V z 



exp 



\_-p_ 
P + l 



[j + l- a) 



(vpr (i+1) E j m,+d 



Z=— oo 



p+l 



(J + 1 - «) 



(x/p)' 
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we obtain from (ITT]) 



M(x) ~ MW(x) = (1 - P ) P c - a+1 J\^J2 i3 + \ 



x exp 



x(A + p) l-3p 
■^H + — —rh + l-a) 



8.2 Approximation for x = 0(e 2 ) = O (c 2 ) 

We shall now use the asymptotic solution in the region R, as given by fl93|) . We 



oo 

M(x) ~ M^(y) = -VG (y, fce) -= / exp 



~K(y, z)dz. 



To evaluate ( II 3 II) as e — > 0, we use the Laplace method. From ( 1201) we get 

^ 2 (y,2;) = g = <S> st = ln(p) 
and therefore the main contribution to ( I131I) comes from z = 1 and we obtain 



M^\y) 



Using (TTTdD in CT-flSgj) and (J^-®, we obtain 
<MV,1) = ^£ K(y,l) 



exp 



1 



K(3/,l) 



1-p /-%,!) 



while (1211) gives 

with 
Thus, 



/' A ' " ' ' ' 2.S'(//. 1 ) \: ••// 
*( 2 /,l) = 2 3 /5 , ( 1 / J l) + lii(p), 



C= v / 2(A-/i)-(A + / i)ln(p). 



M {2 \y) 



p — A 1 — p 



exp 



;-2Cv^ + lnp) 
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9 Summary and discussion 

In most of the domain Q = {(y, z) : y,z > 0} , the asymptotic expansion of Fk(x) = G(y, z) 
is given by 



G(y, z) ~ v/iexp 



e 



K{y,z) mR 



c 



or 



G(oo, z) — G(y, z) ~ — y^exp 



~^{y,z) 



K(y, z) in R. 



(136) 
(137) 



If we consider the continuous part of the density, given by 

dG 



f k (x) = F^x)=e 2 — (y,z), x>0 



dy 



the transition between R and R disappears, and we have 



f k {x) ~ e^ y (y,z)exp 



K(y, z) = e 2 exp 



sK{s,t), 



'138) 



everywhere in the interior of D. Note that K(y,z) becomes infinite along y = Yq(z) (i.e., 
s — 0), but the product ^ y (y, z)K(y, z) remains finite. 

The asymptotic expansion of the boundary probabilities F k (0), k < [_cj can be obtained 
by setting y = in (j!37p . This expression can be used to estimate the difference 



F fc (oo) - F fc (0) = Pr X(oo) > 0, Z(oo) = k 



which is exponentially small for e — > 0. Also, for a fixed z G [0, 1), fk{x) is maximal at x = 
(see Figure [3]). In other words, if k < c, the buffer will most likely be empty. 

For a fixed z > 1, fk(x) is peeked along the curve y = Y (z) (see Figure Hj). To see this 
better, we can use ( 1271) . fl36|) and fl37j) in f 1 1 3 8 [) . obtaining 

/*(*) ~ ^^^= v 7 ^) ex P | - [ zln ip) ~ x{z) (y - Y f] I , z > l 



with 



or equivalently 



3 (/i - A) ; 



(/i + A) (z-1) 



3 ' 



cxp 



a; - c 2 Yn 



^ > 1. 
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Figure 3: A plot of e^ y (y, z) exp [\^{y,z]] K(y,z), with e = 0.1, A = 0.3145, \i = 0.8473 
and z = 0.5. 



31 




32 



4e-09 






3e-09 


/ □ n □ \ 
/ □ □ \ 

/ ° n \ 




2e-09 


/ a \ 
/ □ \ 

/a n \ 

/ U u \ 
/ □ \ 

/ n \ 
/ □ \ 




1e-09 - 


/ U a \ 
/ a + * + + + * + \ 
/ □ + + + n \ 
/ □ + + \ 

/ D ° \ 
! S S o o ° ° ° ° o s^i-g 







0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

z 


2 



Figure 5: A plot of eitf tf (y, z) exp K(j/,3), with e = 0.1, A = 0.3145,// = 0.8473 

for y = 0.5 (solid line), y = 0.6 (□), y = 0.8 (+++) and y — 1 (ooo). 

This means that given k > c active sources, the most likely value of the buffer will be 



For a fixed x > 0, /^(rr) achieves its maximum around z = 1 (see Figure [5]). To find an 
expression for fk{x) when z is close to 1, we use (11331) and obtain, for fixed y > 0, 



r / x 3l-p/C_3 jl 



lnp-2C^ 



or 



fk(x) ~ ^ (1 - p) P C \J^ X * ex P 



C(^-i) 2 

2^(/i-A) 

C(fc-c) 2 
(/i — A) 



Below we summarize the various boundary, corner and transition layer corrections to the 
results in (j93|) and (|92|): 
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1. k = l + c-a, x = 0(l) : 



Fp\x) = (l-p) 



' p — X 



.c-a+i 



n + x 

X + p 



Br 



cxp [A(0)] dO, 



where J.(-) denotes the Bessel function, T (•) the Gamma function, Br is a vertical 
contour in the complex plane with Re(s) > and 

A + p 



a = c-[c|e(0,l), P=^<1, A(0) = y-(^-a)ln 



\/P 



e 



2. y-FoW = 0(v^), Kz: 



F fc (:r)~(l-p)p fe - 



1 + erf ( 



with 

V(y,z) 
3. fc = 0(1) 



^/e V 2 // + A 



F fc (3) (y)-F fc (oo)~ v^exp 



2(/i-A)' 



:*(y,0) 



X 



a ' y) ^(y) + P fc r fc (y)' 



pi — X 



27r3 Q (y)S(y,0) 



(1-P), 



p£(?/) - 



A 



y-KO, 



0) = 2yS(y, 0) + In < 0, J (y) = 2 

5( 2/ ,0) = (A + /x)-^(y)-Ar 1 (y)<0, 

(i - r 1 ) p - a - o - (p + 1) in (o = I* [(i - r 1 ) p + (i - o] 2 y- 

A. y = eu, u — 0(1), 1 < z 

(X + p)u 



F k (x) ~ e (1 - p) — r 

V ^ V ^ V 1+P27T2-1 

f In (p) Xe 2 eu 

x exp ^ — — H In 

[ e e 



z-1 



(z-iy 



+ a In 



(A + /i)w, 

x exp < — In 

z — 1 



(p+l)(^-l) 



+ 



+ 1 - a 

(A + /i) u 

2- 1 
2Aw 



z-1 
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